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ABSTRACT 


% 


Within  the  stratus-topped  boundary  layer  there  axe  many  physical  processes  involved: 
longwave  radiation  cooling,  entrainment,  latent  heating,  surface  heating,  solar  heating, 
drizzling,  etc.  How  all  processes  combine  to  maintain  the  turbulence  within  the  stratus- 
topped  boundary  layer  remains  an  unsolved  problem.  The  large-eddy  simulation  tech¬ 
nique  is  used  to  examine  the  first  four  physical  processes  mentioned  above.  First,  the  con¬ 
tribution  of  each  physical  process  to  the  thermodynamic  differences  between  the  updraft 
and  downdraft  branches  of  turbulent  circulations  is  examined  through  a  conditional  sam¬ 
pling.  Second,  these  mean  thermodynamic  differences  are  shown  to  express  well  the  verti¬ 
cal  distributions  of  heat  and  moisture  fluxes  within  stratus-topped  boundary  layers. 

These  provide  a  method  to  validate  the  process  paxtitioning  technique.  (This  tech¬ 
nique  assumes  that  the  net  flux  profile  can  be  partitioned  into  different  component-flux 
.  profiles  according  to  physical  processes,  amd  that  each  partitioned  component-flux  is  lin¬ 
ear  in  height.)  In  this  paper,  the  heat  and  moisture  fluxes  are  process-partitioned,  and 
,  ’each  component-flux  is  found  to  contribute  to  the  net  flux  in  a  way  that  is  consistent  with 
its  corresponding  process  contribution  to  the  mean  thermodynamic  differences  between 
updrafts  and  downdrafts.  Also,  the  net  flux  obtained  by  summing  all  component-fluxes 
agrees  very  well  with  that  obtained  directly  from  the  large-eddy  simiilations. 
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1.  Introduction 


The  nocturnal,  stratus-topped  boundary  layer  has  many  important  physical  processes: 
longwave  radiation,  entrainment,  surface  heating,  latent  heating,  drizzling,  etc.,  all  work¬ 
ing  interactively  and  simiiltaneously  in  maintaining  the  turbulent  circulations  that  trans¬ 
port  most  of  the  momentum,  heat,  and  moisture.  How  these  processes  combine  to  de¬ 
termine  the  turbulent  tramsports  is  poorly  understood,  however,  mainly  for  the  follow¬ 
ing  reasons:  (1)  Most  of  these  processes  occur  (and  interact)  near  the  cloud  top  where 
the  turbulent  boundary  layer  and  nonturbulent  free  atmosphere  meet.  (2)  The  strengths 
and  the  roles  of  these  processes  and  their  interactions  vary  in  location  along  the  cloud  top 
(NichoUs,  1989;  Moeng  and  Schumann,  1991).  For  example,  more  dry  inversion  air  is  en¬ 
trained  (i.e.,  engudfed)  into  downdrafts  of  major  circulations,  resulting  in  more  evapora¬ 
tive  cooling  there  than  in  updrafts.  This  kind  of  information  is  lost  when  examining  typi¬ 
cal  statistics,  i.e.,  quantities  that  are  defined  in  terms  of  area  averages. 

Furthermore,  careful  and  detailed  observational  studies  of  marine  stratus  ue  rare. 

One  reason  is  that  typical  marine  stratus  is  not  normally  fotmd  near  shore,  since  conti¬ 
nental  effects  can  alter  the  boundary-layer  structure  and  disrupt  cloud  persistence.  As  a 
result,  most  observational  studies  have  relied  heavily  on  long-range  aircraft,  which  can¬ 
not  provide  much  information  on  time  evolution  or  the  three-dimensional  flow  patterns. 
This,  and  the  fact  that  observational  data  are  affected  by  large-scale  horizontal  and  diur¬ 
nal  variability,  and  by  decoupling  of  cloud  layers,  makes  it  difficult  to  study  [or  isolate]  the 
physical  processes,  their  interactions,  and  their  roles. 

In  this  study,  we  use  the  large-eddy  simulation  (LES)  approach  (Moeng,  1984, 1986) 
for  studying  turbulent  flows,  which  was  first  pursued  by  Deardorff  in  the  early  1970s,  and 
has  resulted  in  major  advances  in  our  imderstanding  of  the  cleu  convective  boundziry 
layer  (e.g.,  Deardorff,  1974a,  b;  Moeng  and  Wyngaard,  1984, 1986, 1989).  The  gmdelines 
for  its  application  to  planetary  boundary-layer  (PBL)  research  were  reviewed  by  Wyn¬ 
gaard  (1984).  Deardorff  (1980a)  and  Moeng  (1986,  1987)  applied  the  LES  approach  to 
stratus-topped  bormdary  layers.  These  simulations  not  only  provide  three-dimensional, 
time-dependent,  stochastic  fields  of  velocity,  temperature,  and  moisture  for  detailed  anal¬ 
ysis  but  also  can  be  “controlled”  so  we  can  systematically  examine  each  bsisic  process  in¬ 


volved. 
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In  Section  2,  we  describe  the  large-eddy  simulated  cases.  In  Section  3,  we  review  the 
convective  mass  fltix  approach  to  boundary  layer  modeling,  and  examine  the  relationship 
between  the  turbulent  fluxes  and  the  thermodynamic  differences  between  updraft  and 
downdraft  branches  of  the  turbtdent  circulations.  We  then  study  the  roles  of  several  im¬ 
portant  physical  processes  in  determining  these  updraft  and  downdraft  properties  (Section 
4),  and  in  determining  the  vertical  distributions  of  the  total  thermodynamic  energy  and 
moisture  fluxes  (Section  5).  We  summarize  and  conclude  this  study  in  Section  6. 

2.  The  large-eddy  simulations 

Three  large-eddy  simulations  of  idealized  nocturnal  stratus-topped  boundary  layers 
have  been  carried  out  to  study  the  relevant  physical  processes.  These  simulations  used  80 
X  80  X  80  grids  covering  a  numerical  domain  of  5  km  x  5  km  x  1  km.  Therefore,  they 
resolve  turbulent  eddies  of  sizes  down  to  125  m  x  125  m  x  25  m,  twice  the  grid  spac¬ 
ing.  Moeng  and  Schumann  (1991)  have  analyzed  these  LES-data,  but  they  studied  just 
the  composite  coherent  structures.  In  this  paper,  we  further  analyze  these  simulations  to 
determine  the  relevant  physical  processes  for  generating  and  maintaining  the  turbulent 
transports.  We  discuss  mainly  (1)  cloud-top  longwave  radiation,  (2)  entrainment,  (3)  sur¬ 
face  heating,  and  (4)  latent  heating,  which  we  believe  are  the  most  important  mechanisms 
for  maintaining  or  dissipating  the  nocturnal  marine  stratus.  (Drizzling  may  also  be  impor¬ 
tant  in  dissipating  the  cloud,  but  we  will  not  consider  that  process  here.) 

The  three  simulations  of  idealized  stratus-topped  PBLs  ue  summarized  in  Table  1. 
The  ENT/RAD  simulation  includes  mainly  the  entrainment  and  longwave  radiation  pro¬ 
cesses.  In  this  case,  we  minimized  the  surface  heat  flux.  (Since  we  prescribed  the  sea  sur¬ 
face  temperature  in  the  simulations,  the  surface  heat  flux  caxmot  be  exactly  zero,  as  we 
wish.)  We  also  excluded  latent  heating  by  using  the  '*dust  cloud”  concept,  introduced  by 
Lilly  (1968).  In  a  dust  cloud,  “liquid  water”  can  be  “condensed”  from  or  “evaporated”  to 
water  vapor,  but  no  latent  heating/cooling  occurs.  The  ENT/RAD/SFC  simulation  in¬ 
cludes  entrainment,  longwave  radiation,  and  s\irface  heating;  latent  heating  is  again  ex¬ 
cluded  by  considering  a  “dust  cloud.”  The  ENT/RAD/SFC/EVP  simulation  includes  all 
four  processes,  and  thus  represents  a  realistic  nocturnal,  stratus-topped  boundary  layer. 


4 


This  simuiation  has  been  used  by  Schumann  and  Moeng  (1991a,  b)  to  study  plume  fluxes 

and  budgets,  and  by  Rajidall  et  al.  (1992)  to  test  a  mass  flux  model. 

We  can  then  study  the  effects  of  evaporation  by  comparing  the  ENT/RAD/SFC 
and  ENT/RAD/SFC/EVP  cases,  and  the  effects  of  surface  heating  by  comparing  the 
ENT/RAD  and  ENT/RAD/SFC  cases.  It  is  more  difficult  to  separate  the  radiative  cool¬ 
ing  and  entreiinment  processes,  since  they  always  coexist  in  even  the  most  idealized,  pro¬ 
totype  stratus  cases.  Therefore,  we  include  in  this  study  an  LES  of  the  clear  convective 
boundary  layer  (CBL)  that  is  driven  mainly  by  surface  heating.  Since  only  surface  flux 
and  entrainment  processes  are  involved,  we  label  it  ENT/SFC  (also  summarized  in  Table 
1).  The  ENT/SFC  case  has  been  studied  by  Moeng  and  Wyngaard  (1988, 1989),  Schu¬ 
mann  and  Moeng  (1991a,  b),  Moeng  and  Schumann  (1991),  and  Wyngaard  zind  Moeng 
(1992). 

The  statistics  given  in  Table  1  are:  the  average  boundary-layer  top  height,  Zt;  the 
average  cloud-base  height,  zb]  the  friction  velocity  u«;  the  surface  buoyancy  flux  WSZq] 
the  convective  velocity  scale  it;,  =  Wc»  =  for  the  stratus  cases,  and  w»  = 

{j^ZivU^oY^^  for  the  clear  boundary  layer  case;  the  temperature  scale,  T,  2  B»/we»  * 
for  the  stratus  cases,  and  T«  =  vjs^Q/{cpW^)  for  the  clear  case;  and  the  moisture  scaling, 
Q,  =  M^/wen  for  the  stratus  cases  only.  Here  J?,  =  ^  vJaZdzf^Zi  —  zb)  and  Af,  = 

//*  Ti;(gv  +  <ll)dz/{zi  —  Zb)  are  the  vertically  averaged  buoyancy  and  moisture  fluxes  within 
the  cloud  layer;  g/T^  is  the  buoyancy  coefficient;  Cp  is  the  specific  heat  of  air;  s,  is  the  vir¬ 
tual  dry  static  energy;  hi  is  the  liquid  water  static  energy,  a  conserved  quantity  under  both 
moist  and  dry  adiabatic  processes  (Moeng,  1986);  g,  is  the  water  vapor  mixing  ratio;  qi 
is  the  liquid  water  mixing  ratio;  and  w  is  the  vertical  velocity.  The  s,  and  hi  fields  have 
the  dimensions  of  temperature  after  we  divide  them  by  Cp.  In  the  ENT/SFC  simulation, 
moisture  was  not  included;  there,  the  predicted  thermal  variable  is  s,,  which  is  conserved 
under  dry  adiabatic  processes.  (The  statistics  we  show  here  for  the  ENT/RAD/SFC/EVP 
case  are  from  a  later  time  period  than  the  dust  cloud  cases,  so  the  cloud  top  is  higher.  We 
did  not  save  the  LES  data  of  the  earlier  simulation  period  for  the  ENT/RAD/SFC/EVP 
case.) 

The  total  (solid  curves)  and  resolved  (dashed)  buoyancy  fluxes  for  these  four 
large-eddy  simulations  are  given  in  Fig.  1.  The  upper-left  one  is  ENT/RAD  (dust 
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cloud  case),  the  lower-left  is  ENT/RAD/SFC  (dust  cloud  case),  the  upper-right  one  is 
ENT/RAD/SFC/EVP  (wet  cloud  case),  and  the  lower-right  is  ENT/SFC  (clear  PEL 
case).  The  dust  or  wet  cloud  layers  are  marked.  The  fluxes  are  shown  in  units  of  W  m“^, 
and  the  height  is  normalized  by  z,-.  The  statistics  are  averages  over  horizontal  surfaces  and 
over  a  tiu-e  period  of  about  2  ~  3  large-eddy  turnover  times  (six  LES  time  records).  (The 
local  maximum  above  the  boundary-layer  top  is  likely  due  to  numerical  spurious,  associ¬ 
ated  with  diiflculty  of  the  advection  scheme  to  handle  the  shzirp  temperature  and  mois¬ 
ture  gradients;  the  longer  the  simulation  period  is,  the  larger  the  spurious  maximum  be¬ 
comes.)  Simulations  ENT/RAD/SFC  and  ENT/RAD/SFC/EVP  have  their  surface  buoy¬ 
ancy  fluxes  (generated  by  surface  heating)  comparable  to  their  near  cloud-top  buoyancy 
fluxes  (generated  by  cloud- top  radiative  cooling.)  Note  that  the  ENT/SFC  case  has  a  very 
strong  surface  buoyancy  flux:  240  W  m“*. 

The  ENT/RAD/SFC/EVP  (wet  cloud)  case  has  a  smaller  surface  buoyancy  flux 
(and  a  smaller  surface  moisture  flux  as  well)  than  the  ENT/RAD/SFC  (dust  cloud) 
case  because  of  the  following  reason.  The  sea  surface  temperatures  are  the  same  in  both 
cases.  Since  the  LES  field  of  ENT/RAD/SFC/EVP  (wet  cloud  case)  we  analyze  is  from' 
a  later  time  period  than  other  two  stratus  cases,  it  has  a  warmer  and  wetter  PEL  than 
ENT/RAD/SFC  (dust  cloud  case),  thus  smaller  surface  fluxes.  Although  the  LES  simu¬ 
lations  have  reached  a  quasi-steady  state,  the  mean  temperature  and  moisture  within  the 
mixed  layer  Me  still  changing  in  time.  Therefore,  the  difierences  between  ENT/RAD/SFC 
and  ENT/RAD/SFC/EVP  are  not  totally  due  to  the  latent  heating  effect.  We  will  note 
the  differences  that  are  resulted  from  this  different-time-period-averaging,  wherever  it  ap¬ 
plies. 

3.  The  convective  mass  fluxes 

Convective  mass  flux  modeling,  first  applied  by  Arakawa  (1969)  and  Betts  (1976)  in 
cumulus  parameterization,  is  based  on  the  idealization  that  convective  circulations  consist 
of  two  branches,  updrafts  and  downdrafts,  and  that  the  turbulent  fluxes  associated  with 
these  circulations  can  be  expressed  by 


^  =  Mcifu  -  /<<)> 


(1) 
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where  Me  =  a(l  —  a){wu.  —  Wd)  is  the  convective  mass  flux,  a  is  the  area  covered  by  up¬ 
drafts,  TOu  and  vjd  axe  the  mean  vertical  velocities  of  the  updraft  and  downdraft  branches, 
respectively,  «uid  and  /„  are  the  mean  values  of  any  scalars  within  the  updraft  and 
downdraft  (e.g.,  R£uadall,  1987),  This  approach  was  recently  applied  to  PBL  parameter¬ 
ization  by  Chatfield  and  Brost  (1987),  Randall  (1987),  Wang  md  Albrecht  (1990),  and 
Rand£dl  et  al.  (1992).  One  of  the  closures  in  this  type  of  modeling  is  determining  JV/e. 

In  developing  techniques  for  measuring  the  turbulent  fluxes,  Businger  and  Oncley 
(1990)  proposed  a  “relaxed  eddy  accumulation”  method,  in  which  the  vertiezd  turbulent 
flux  of  a  scalar  species  /  in  the  PBL  is  taken  to  be 

to/  =  6<r«,(/tt  -  fd),  (2) 

where  b  is  a  coefficient  of  proportionality  and  is  the  root-mezm-square  of  the  vertical 
velocity  field.  From  data  collected  in  the  surface  layer,  Businger  and  Oncley  found  that  b 
is  about  0.6,  independent  of  stability. 

Wyngaard  and  Moeng  (1992),  using  the  joint  probability  density  distribution  of  to  ahd 
a  scalar  /  generated  from  LES  data,  computed  the  Businger-Oncley  parameter  b  in  (2). 
They  foxmd  6  ~  0.6  for  a  bottom-up  scalar  (i.e.,  a  scalar  with  zero  entrainment  flux)  and 
0.47  for  a  top-down  scalar  (i.e.,  a  scalar  with  zero  surface  flux)  throughout  the  whole 
CBL.  They  also  showed  analytically  that  for  jointly  Gaussian  turbulence,  h  —  \/^/4  ~ 
0.627. 

Moeng  and  Schumann  (1991)  show  that  turbulence  driven  by  cloud-top  radiative  cool¬ 
ing  has  a  less  skewed  vertical  velocity  field  than  that  driven  by  surface  heating.  We  there¬ 
fore  expect  that  the  joint  probability  density  of  w  and  any  scalar  is  more  jointly  Gaussian 
(i.e.,  b  0.6)  for  the  stratus-topped  PBL  than  for  the  CBL.  We  generate,  firom  the  LES 
results,  the  joint  density  of  w  and  hi,  and  of  w  and  gxt  where  gr  is  the  total  moisture  mix¬ 
ing  ratio.  We  found  that  both  joint  density  distributions  are  indeed  approximately  joint 
Gaussian,  and  b  ~  0.6  throughout  most  of  the  PBL.  Instead  of  showing  the  b  distribution, 
we  compare  whi  (or  wJ^  for  the  clear  case)  obtained  by  two  methods:  (a)  directly  corre¬ 
lating  the  resolved-scale  w  and  hi  (or  for  the  clear  case)  fields  of  the  LES  data  and  (b) 
using  Eq.  (2)  with  6  =  0.6.  The  result  is  shown  in  Fig.  2;  the  agreement  is  remarkably 
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good  for  all  ca^es.  (We  exauiine  the  resolvable  field  only  since  we  cannot  split  the  subgrid 
scales  into  updrafts  and  downdrafts.) 

Figure  3  compares  obtained  (a)  directly  from  the  LES  data  and  (b)  by  using  Eq. 
(2)  with  b  =  0.6.  Again,  the  parameterized  moisture  flux  is  remarkably  good.  No  moisture 
data  are  available  for  the  CBL  case. 

As  shown  by  Randall  et  al.  (1992),  with  the  “top-hat”  model  the  root-mean-square 
vertical  velocity  is  given  by  =  ^0(1  —  a)(ti;u  —  Combining  this  with  Me  = 

a(l  —  a){w^  —  Wd)  gives  b  =  \Ja(l  —  a).  With  this  formula,  the  largest  possible  value  of  b  is 
0.5.  This  implies  that  if  we  use  the  top-hat  model  and  require  consistency  between  (1)  and 
(2),  the  difference  /„  —  /</  between  the  updrzfft  and  downdraft  thermodynamic  variables 
will  have  to  be  exaggerated  by  the  top-hat  model  in  order  to  produce  the  correct  fluxes. 
Further  discussion  is  given  by  Randall  et  al.  (1992). 

The  turbulent  transports  of  heat  and  moisture  are  thus  well  expressed  by  the  different 
properties  of  the  updraft  and  downdraft  branches.  To  see  how  the  four  processes  under 
study  here  contribute  to  these  turbulent  transports,  it  is  necessary  to  examine  how  they 
determine  the  properties  within  the  updraft  and  downdraft  branches  separately. 

4.  Mean  properties  within  the  updrafts  and  downdrafts 
a.  Analyau  methods 

The  local  cloud-top  height,  ZT(a;,y),  varies  in  time  and  space.  Using  conventional  hor¬ 
izontal  averaging,  the  statistics  near  the  cloud-top  level  include  both  above-cloud  and  in¬ 
cloud  properties.  It  is  difficult  to  retrieve  just  the  in-cloud  properties  from  such  an  aver¬ 
aging.  We  therefore  perform  a  new  averaging  process.  We  scale  all  variables  at  each  LES 
grid  colimm  with  the  local  cloud-top  height.  That  is,  we  first  interpolate  the  variables 
at  each  (x,  y)  grid  column  into  a  new  vertical  coordinate  that  puts  all  the  local  cloud-top 
heights  at  the  same  level  in  the  new  vertical  grids.  We  then  perform  “horizontal”  averag¬ 
ing  over  the  new  vertical  coordinate,  which  consists  of  “hilly”  surfaces  that  are  parallel  to 
the  “hilly”  cloud-top  surface.  This  type  of  averaging  was  mentioned  by  Lilly  and  Schubert 
(1980),  and  used  by  Nieuwstadt  and  Businger  (1984).  The  new  “horizontally”  averaged 
statistics  will  be  slightly  different,  especially  near  the  PBL  top,  from  conventional  statis¬ 
tics  calcTilated  by  averaging  over  fiat  horizontal  surfaces. 
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It  is  harder  to  perform  this  type  of  averaging  in  the  clear  convective  case,  since  there 

is  no  cloud  top  to  serve  as  a  reference.  Instead,  we  use  the  temperature  profiles  at  each 
grid  column;  the  local  PBL  top,  ZT[x,y),  is  defined  as  the  height  where  the  temperature 
increases  by  at  least  0.16  K  across  one  grid  interval  ( Az  ~  20  m).  The  average  of  zt{x,  y) 
defined  this  way  equals  approximately  Zi  (the  height  of  the  minimum  buoyimcy  flux). 

We  use  this  averaging  method  only  in  this  section  to  show  the  up  draft /downdraft 
properties.  Statistics  represented  by  this  new  “horizontal”  averaging  will  be  presented 
in  figures  that  have  their  height  coordinate  as  z/zx{x,y).  In  other  sections  of  this  paper, 
statistics  are  represented  by  conventional  horizontal  averaging,  and  the  figiires  use  the 
z/zi  height  coordinate. 

We  also  perform  conditional  sampling  to  separate  the  updraft  and  downdraft  branches 
of  the  circulations.  The  vertical  velocity  fluctuations  are  used  as  an  indicator:  lu  >  0  for 
the  updraft  branch  and  w  <  0  for  the  downdraft  branch.  This  conditional  sampling  tech¬ 
nique  has  been  used  by  many  others  (e.g.,  Greenhut  and  Khalsa,  1987;  Young,  1988;  Schu¬ 
mann  and  Moeng,  1991a). 

Some  research  (Khalsa,  S.J.S.,  personal  communication,  1991)  suggests  that  positive, 
vertical  velocity  fluctuations  near  the  PBL  top  do  not  necessarily  coincide  with  the  up¬ 
drafts  of  major  circ\ilations.  Therefore,  using  local  vertical  velocity  as  *ux  indicator  for 
splitting  the  updraft  and  downdraft  branches  of  the  dominant  circulations  may  be  inap¬ 
propriate.  Thus,  we  also  perform  a  diflerent  type  of  conditional  sampling  that  splits  the 
updraft  and  downdraft  branches  according  to  the  vertical  velocity  fluctuations  in  the  mid¬ 
dle  of  the  PBL.  Even  though  the  results  (not  shown)  differ  quantitatively  from  those  using 
the  local  vertical  velocity  fleld  as  the  indicator,  the  resxilts  presented  below  are  applicable 
to  this  type  of  sampling  as  well.  In  the  following,  our  definition  of  the  updraft  branch  is  a 
collection  of  all  the  ascending  elements,  and  that  of  the  downdraft  branch  a  collection  of 
all  the  descending  elements. 

b.  Results  and  discussion 

The  solid  and  dashed  curves  in  Figs.  4,  5,  6,  and  7  show  the  vertical  distributions  of 
the  mean  properties  within  the  updraft  branch  and  downdraft  branch,  respectively.  The 
four  panels  are  for  the  different  LES  cases.  These  figures  show,  in  order,  the  vertical  ve- 
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locity,  virtual  temperature,  total  moisture,  and  liquid  water  mixing  ratio  that  are  "hori¬ 
zontally”  averaged  (over  the  “hilly”  surface)  within  updrafts  and  within  downdrafts,  re¬ 
spectively.  The  “horizontal”  means  that  are  subtracted  out  in  these  figures  have  similar 
profiles  to  the  conventional  horizontal-mean  profiles;  the  latter  for  ENT/RAD/SFC/EVP 
(wet  cloud)  and  ENT/SFC  (clear  PEL)  cases  are  reported  in  Schumann  amd  Moeng 
(1991a). 

The  vertical  velocity  differences  between  updrafts  and  downdrafts  for  all  three  stratus 
cloud  cases  are  within  30%;  the  addition  of  surface  heating  to  cloud-top  radiative  cooling 
increases  the  circulations  (measured  as  —  Wd)  by  about  30%  in  the  lower  PEL.  The  cir- 
ctdations  in  the  ENT/SFC  (clear  PEL)  case  are  particularly  strong,  since  the  input  sur¬ 
face  buoyancy  flux  is  particularly  large  (Fig.  1).  Updrafts  are  about  twice  as  strong  as 
downdrafts  in  the  ENT/SFC  case;  i.e.,  the  CEL  has  a  positively  skewed  vertical  velocity 
field,  as  discussed  in  Moeng  and  Schumann  (1991). 

Figure  5  shows  the  mean  virtual  temperature  differences  between  the  updrafts  and 
downdrafts.  Within  the  upper  5%  of  the  dust  cloud-topped  PEL  (i.e.,  in  the  ENT/RAD 
and  ENT/RAD/SFC  cases  where  the  evaporation/ condensation  effect  is  neglected),  the 
updraft  branch  is  colder  than  the  downdraft  branch.  [Note  that  within  this  5%  layer  the 
air  is  totally  cloudy  since  we  scale  the  height  with  r7(x,y).]  In  the  upper  part  of  the  up¬ 
draft,  there  is  less  entrained  inversion  air  and  also  there  is  adiabatic  lifting;  both  con¬ 
tribute  to  make  the  updrafts  colder  than  the  downdrafts  in  this  uppermost  layer.  When 
evaporative  cooling  is  included,  however,  as  in  the  ENT/RAD/SFC/EVP  (wet  cloud) 
case,  the  downdrafts  become  colder  than  updrafts  everywhere.  Apparently,  cloud-top 
evaporation  cools  downdrafts  more  than  updrafts,  for  more  dry  air  is  engulfed  into  down- 
drafts.  This  latent-heating  contribution  to  the  up  draft /downdraft  differences  offsets  the 
net  effect  from  entrainment  and  adiabatic  lifting,  resiilting  in  colder  downdrafts.  Later,  we 
will  use  the  liquid-water  field,  given  below,  to  estimate  this  contribution. 

Within  the  mixed  layer,  updrafts  are  warmer  than  downdrafts,  so  the  turbulent  cir¬ 
culations  are  buoyantly  driven.  The  s\irface  heating  in  the  ENT/RAD/SFC  (dust  cloud) 
case  results  in  much  larger  temperature  variations  in  the  lower  half  of  the  PEL  than  in 
the  ENT/RAD  (dust  cloud)  case;  consequently,  the  heat  fiux  is  much  larger  there,  as  will 
be  shown  in  Section  5.  The  latent  heating  causes  smaller  temperature  variations  near 
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the  cloud  base  in  ENT/RAD/SFC/  EVP  (wet  cloud)  than  in  the  ENT/RAD/SFC  \^dust 
cloud)  case.  This  seems  to  imply  that  condensation  heats  downdrafts  more  than  updrafts, 
which  is  puzzling.  (The  difference  in  the  lower  PEL  between  the  ENT/RAD/SFC  and 
ENT/RAD/SFC/EVP  cases  is  mainly  due  to  the  fact  that  their  statistics  are  taken  from  a 
different  time  period.) 

The  updraft/downdraft  temperature  differences  are  more  height-dependent  in  the 
ENT/SFC  (clear  PEL)  case  than  in  the  stratus  cases  because  of  entrainment.  En¬ 
trainment  brings  more  warm  air  into  downdrjffts  than  updrafts,  and  reduces  the  up¬ 
draft/downdraft  temperature  differences  from  surface  heating.  This  effect  is  greatest  at 
the  PEL  top  and  decreases  downward,  while  the  surface  heating  effect  has  a  maximum  at 
the  surface  and  decreases  upward.  For  the  stratus  cases,  the  reduction  of  (r„)u  —  iTv)d  due 
to  entrainment  works  concurrently  with  the  enhancement  of  (r„)u  —  {Tv)d  due  to  cloud-top 
radiative  cooling,  thus  resulting  in  less  height  dependence  of  (T«)u  —  (T'«)a. 

Figure  6  shows  the  mean  moisture  differences  between  the  updrafts  and  downdrafts. 
The  updrafts  are  always  wetter  than  the  downdrafts  throughout  the  whole  PEL,  as  a 
result  of  both  surface  and  entrainment  processes.  (The  radiation  and  evaporation  pro-  * 
cesses  do  not  contribute  to  the  total  mixing  ratio.)  Even  though  the  surface  buoyancy 
flux  is  negligibly  small  in  the  ENT/RAD  (dust  cloud)  case,  there  is  a  significant  sur¬ 
face  moisture  flux  because  the  saturated  mixing  ratio  at  the  ocean  surface  is  much  Isirger 
th8m  the  mixing  ratio  within  the  PEL.  Therefore,  the  moisture  differences  between  up¬ 
drafts  and  downdrafts  in  the  ENT/RAD  (dust  cloud)  case  are  due  to  both  entrainment¬ 
drying  and  surface-moistening  processes.  Including  the  surface-heating  process  in  the 
ENT/RAD/SFC  (dust  cloud)  case  results  in  moisture  differences  almost  twice  as  large  as 
in  the  ENT/RAD  (dust  cloud)  case.  Latent  heating  seems  to  have  little  effect.  (The  dif¬ 
ference  in  the  lower  PEL  between  the  ENT/RAD/SFC  and  ENT/RAD/SFC/EVP  cases 
is  mainly  due  to  the  fact  that  their  statistics  are  taken  from  a  different  time  period.) 

The  liquid-water  mixing  ratio  (Fig.  7)  is  larger  within  the  upper  10%  of  the  updrafts 
than  downdrafts.  This  is  due  to  a  larger  total  mixing  ratio  in  the  updrafts.  The  differ¬ 
ences  in  qi  can  contribute  to  the  evaporative  cooling  only  in  the  ENT/RAD/SFC/EVP 
(wet  cloud)  case,  where  the  latent  heating  effect  is  included.  We  can  estimate  the  tem¬ 
perature  differences  between  the  updrafts  and  downdrafts  due  to  evaporative  cooling  as 
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follows:  We  assume  that  the  updraft /downdraft  temperature  differences  due  to  evapora¬ 
tive  cooling  are  i[(gj)u  -  {(li)d]lcp,  where  {qi)^  -  {qi)^  ~  7  x  10“^  g/Kg  at  the  cloud-top 
level  (from  Fig.  7c).  Thus,  the  evaporation  contribution  alone  gives  ~  0.17  K. 

This  is  more  than  enough  to  bring  the  downdradfs  at  the  cloud  top  in  the  dust  cloud  cases 
(which  are  only  about  0.06  K  warmer  in  ENT/RAD/SFC)  to  nearly  the  saune  tempera¬ 
ture  as  updrafts  in  the  real  cloud  case.  Since  the  liquid-water  differences  vary  between  the 
ENT/RAD/SFC/EVP  (wet  cloud)  aind  ENT/RAD/SFC  (dust  cloud)  cases,  it  is  difficult 
to  exactly  figure  out  the  evaporation  effect.  We  adso  note  that  within  the  mixed  layer,  the 
liquid- water  differences  between  updrafts  and  downdrafts  inexplicably  change  from  case  to 
case. 

The  solid  and  dashed  curves  in  Fig.  8  show  the  mean  longwave  radiative  fluxes  within 
the  updrafts  and  downdrafts,  respectively.  We  compute  (Frad)u  (.frtt<i)<<  by  vertically 
integrating  the  radiative  cooling  rate  at  each  (x,y)  column  and  averaging  over  updrafts 
and  downdrafts,  respectively.  (We  have  saved  the  longwave  radiative  cooling  rate  at  each 
grid  point  in  the  LES  records.)  Well  within  the  cloudy  layer,  we  set  the  radiative  fluxes  to 
zero.  We  multiply  the  flux  by  Cp  to  obtain  \inits  of  W  m~*.  The  mean  longwave  radiative 
fluxes  within  the  updrafts  and  downdrafts  differ  slightly;  the  flux  jump  (i.e.,  the  slope) 
near  the  cloud  top  is  slightly  larger  in  the  updraft  branch  than  in  the  downdraft  brsmch. 

An  interesting  feature  is  that  about  10%  of  the  flux  jump  is  above  the  cloudy  layer 
[i.e.,  above  zj(z,y)],  not  too  different  &om  that  observed  by  Slingo  et  al.  (1982)  and  re¬ 
ported  Nieuwstadt  and  Businger  (1984).  As  we  will  see  in  Section  5,  Ftop — the  mean 
longwave  radiative  flux  averaged  along  the  “hilly”  zt(®>  y)  surface — ^is  about  10%  smaller 
than  the  average  Frad  at  the  Zi  level.  We  define  AF  =  Ftop  —  (i^rai<)o>  which  is  the  bulk- 
flux- jump  across  the  whole  trurbulent  layer,  where  (Fra<i)o  is  tli®  radiation  flux  at  the  sur¬ 
face.  Since  only  the  radiative  cooling  (i.e.,  the  flux  jump)  that  is  entirely  within  the  turbu¬ 
lent  layer  can  drive  the  turb\ilent  circulations  directly  (e.g.,  Kahn  and  Businger  1979),  we 
will  use  AF  in  Section  5  to  partition  the  total  thermodynamic  energy  flux  profile. 

The  radiative  cooling  above  the  cloudy  layer  (i.e.,  above  zx(®>y))  may  cool  the  in¬ 
version  air  and  mak^  the  entrainment  process  easier.  In  that  sense,  the  radiative  cooling 
above  ZT(x,y)  may  be  important  in  determining  the  entrainment  rate,  thus  may  indirectly 
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anect  the  turbuleut  transports.  In  this  paper,  we  do  not  intend  to  address  the  problem 

determining  entrainment  rate,  however. 

c.  Summary 

To  conclude  this  section,  we  summarize  the  roles  of  each  physiced  process  in  determin¬ 
ing  the  mean  properties  within  updrafts  2Lnd  downdrafts,  separately. 

(1)  Moeng  and  Schumann  (1991)  show  that  downdrafts  near  the  cloud  top  axe  ini¬ 
tially  forced  as  the  ho:izontal  flow  converges.  As  downdrafts  axe  forced  downward  into 
the  mixed  layer,  more  radiatively  cooled  air  paxcels  tire  incorporated  into  them  than  into 
updrafts.  This  way,  cloud-top  radiative  cooling  contributes  to  colder  downdrafts  than 
updrafts,  thus  bu  /antly  driving  the  whole  turbulent  circulations.  The  about  90%  of  the 
cloud-top  radiative  flux  jump  that  is  within  the  cloudy  layer  can  actuzdly  drive  the  turbu¬ 
lent  circulations  2uid  contribute  directly  to  the  turbulent  fluxes,  as  we  will  show  in  Section 
5. 

(2)  Entrainment  brings  warm  and  dry  inversion  air  into  the  PEL,  mostly  into  the 
downdrafts.  Thus,  entrainment  tends  to  make  the  downdrafts  warmer  than  the  updrafts 
throughout  the  whole  depth  of  the  PEL.  In  the  absence  of  surface  heating,  only  when 
cloud-top  radiative  cooling  dominates  entrainment,  can  the  PEL  be  buoyantly  driven.  En- 
treunment  also  gives  drier  downdrafts  than  updrafts  throughout  the  PEL. 

(3)  Surface  heating  (and  surface  moistening)  seems  to  be  more  efficient,  in  com¬ 
parison  to  the  cloud-top  radiative  cooling  process,  in  producing  warmer  (and  wetter) 
updrafts  than  downdrafts.  From  ENT/RAD  and  ENT/RAD/SFC  (both  dust  cloud 
cases)  we  show  that  the  updrafts/downdrafts  temperature  differences  in  the  lower  half 
of  the  ENT/RAD/SFC  are  more  than  twice  as  large  as  in  the  ENT/RAD  case.  The  up¬ 
draft/downdraft  moisttire  differences  also  increase,  largely  due  to  surface  effects. 

(4)  Evaporative  cooling  is  larger  in  downdrafts  than  updrafts  within  the  upper  5%  of 
the  cloud-top  layer,  since  there  is  more  entrained  dry  air  within  downdrafts.  This  cooling 
overcomes  other  processes,  and  as  a  net  result,  downdrzdts  become  colder  than  updrafts 
also  in  that  layer,  as  shown  in  Fig.  5c.  In  our  simulated  cloud  case,  which  is  stable  with 
respect  to  the  cloud-top  entrainment  instability,  the  latent  heating  contribution  seems  to 
be  confined  to  this  thin  top  layer;  it  seems  to  contribute  very  little  to  the  overall  turbulent 
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circulations,  or  to  the  overall  heat  and  moisture  variations.  When  the  cloud-top  entrain¬ 
ment  instability  condition  is  satisfied,  this  process  may  become  more  significant,  and  a 
strong  evaporation  / entrainment  interaction  may  then  contribute  to  a  breakup  of  the  stra¬ 
tus  deck  (Lilly,  1968;  Randall,  1980a;  DeardorfF,  1980b;  Albrecht,  et  al.  1985;  MacVean 
and  Mason,  1990). 

5.  Contributions  of  processes  to  the  vertical  distributions  of  the  turbulent 
fluxes 

We  have  shown  that  the  turbulent  transports  axe  well  expressed  by  the  differences  be¬ 
tween  the  thermodynamic  properties  of  the  updrafts  and  downdrafts  in  Section  3.  We 
have  also  qualitatively  examined  the  contributions  of  each  physical  process  to  such  dif¬ 
ferences  in  Section  4.  Based  on  that,  we  will  examine  the  linear  “process  partitioning”  of 
the  heat  and  moisture  fluxes.  Process  partitioning,  flrst  proposed  by  Manias  and  Turner 
(1978)  and  labled  by  Randall  (1984),  is  a  parameterization  that  separates  the  roles  of  var¬ 
ious  physical  processes  in  determining  the  flux  profiles.  It  was  used  by  Stage  (uid  Businger 
(1981),  Wilczak  and  Businger  (1983),  and  Hanson  (1987a)  to  partition  the  buoyancy  flux 
within  the  convective  PBL.  It  assumes  that  each  partitioned  fltix  is  lineu  in  height,  and 
the  linear  sum  of  all  partitioned  fluxes  represents  the  total  flux.  Here,  we  will  partition  the 
total  thermodynamic  energy  flux  and  the  total  moisture  flux  since  they  are  linear  in  height 
for  a  homogeneous,  quasi-steady  state  boundary  layer,  and  therefore  are  easier  to  parti¬ 
tion. 

a.  Partitioning  the  thermodynamic  energy  flux 

Figure  9  shows  the  total  thermodynamic  energy  flux,  H  =  pwhi  -f  Frad  (solid  curves) 
and  the  longwave  radiative  flux,  Frad  (dashed  curves)  for  all  three  stratus  cases,  where  p 
the  air  density  is  1  Kg/m®  in  our  simulations.  Here  the  statistics  axe  the  conventional  av¬ 
erages  over  X  —  y  horizontal  planes,  and  over  2  ~  3  large-eddy  turnover  times.  Thus,  F rad 
is  the  horizontally-averaged  radiative  flux.  The  total  thennodynamic  energy  flux  should 
be  linear  in  height  within  the  mixed  layer  (i.e.,  below  ~  0.95zi  level  marked  by  x  in  Fig.  9) 
for  a  homogeneous,  quasi-steady  boundary  layer,  because 
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=  0. 


As  discussed  in  Section  4,  about  90%  of  tbe  radiative  flux  jump  is  witbin  tbe  solid 
cloud  layer  (i.e.,  the  turbulent  layer),  and  is  available  for  driving  the  turbulence  directly. 
Some  of  the  flux  at  Zi  is  actually  above  the  solid  cloud  layer  and  should  be  excluded.  We 
hypothesize  here  that  the  dotted  curves  in  Fig.  9,  which  are  the  linear  curves  that  have 
Ftop  at  Zi  and  (jPra<i)o  at  the  surface,  are  the  J-flux  component  response  to  an  exter¬ 
nal  forcing  of  only  longwave  radiation.  Using  Ftop  (the  average  radiative  flux  along  the 
ZT{x,y)  surface),  instead  of  Fradi^i),  excludes  the  above-cloud  radiation  eflfect  on  turbu¬ 
lent  transports.  For  the  ENT/RAD  (dust  cloud)  case,  where  the  surface  heating  is  negligi¬ 
bly  small,  for  example,  H  is  smaller  than  the  dotted  curve  throughout  the  PEL  because  of 
entrainment.  In  other  words,  entrzunment  reduces  S,  and  thus  whi,  for  a  given  longwave 
radiation  forcing. 

Figure  10  shows  the  vertical  distributions  of  whi  from  the  LES  restilts.  Note  that  they 
are  not  linear  with  height  because  Frai  is  not  linear. 

From  the  above  argument,  we  partition  the  profile  of  B  in  Fig.  11  according  to  the 
various  physical  processes  involved.  We  assume  that  the  contribution  of  each  process  is 
linear  with  height,  as  suggested  by  Stage  and  Businger  (1981)  and  as  implied  from  the 
analysis  of  Section  4.  For  convenience,  we  shift  the  x-coordinates  of  Fig.  9  to  those  of  Fig. 
11  to  eliminate  the  surface  radiation  flux  (which  should  not  aSect  the  turbulent  trans¬ 
ports  anyway).  Note  that  AF  =  Ftop  —  (Fro«i)o*  [AF  is  larger  in  ENT/RAD  than 
in  ENT/RAD/SFC  and  ENT/RAD/SPC/EVP  because  (1)  the  net  longwave  radiative 
fliix  along  the  cloud  top  (i.e.,  Ftop)  is  about  the  same  for  all  three  cases,  and  (2)  the 
ENT/RAD  case  has  a  smaller  sea  surface  temperature  and  thus  a  smaller  (Fra(i)o.] 

The  solid  curves  in  Fig.  11  are  the  B  profiles  in  Fig.  9  after  the  coordinate  shift.  The 
R'-flux  component  response  to  the  cloud-top  radiative  cooling  process  alone  is  represented 
by  the  RAD  curves,  which  are  Ftop  at  Zi  and  zero  at  the  surface.  The  component  re¬ 
sponse  to  the  entrainment  process  alone  is  represented  by  the  ENT  curves,  which  are 

at  Zi  and  zero  at  the  surface.  And,  the  component  response  to  the  surface-heating 
process  alone  is  represented  by  the  SFC  curves,  which  are  (wKi)o  at  the  surface  and  zero 
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at  the  top.  (The  latent  heating  should  not  affect  the  H  profile  since  hi  is  conserved  un¬ 
der  a  moist  adiabatic  process.)  The  sum  (dotted  lines)  of  these  three  curves  gives  approxi¬ 
mately  the  H  profile  computed  from  the  LES  data  (solid  curves). 

Our  process  partitioning  is  basictdly  the  same  as  that  proposed  by  Stage  and  Businger 
(1981),  Hanson  (1987a),  ajid  Stull  (1988,  Fig.  13.20),  and  discussed  by  Randall  (1984). 
However,  since  we  partition  the  H  profile,  not  uJs7  or  whi,  the  vertical  distribution  of  the 
radiative  flux  does  not  affect  the  partitioning.  We  show  that  Frad  is  not  lineax  with  height, 
so  that  whi  (or  the  equivalent  potential  temperature  flux)  is  not  lineeir  either.  The  details 
of  the  Frad  distribution  will  modify  the  whi  profile,  after  we  determine  the  H  profile. 

The  latent  heating  effect  results  in  the  difference  between  whi  and  wJ^,  as  in  the 
ENT/RAD/SFC/EVP  (wet  cloud)  case  shown  in  Figs.  1  &  10.  (These  two  fluxes  axe 
identical  in  the  “dust  cloud”  simulations.)  The  buoyancy  flux  is  of  great  importance, 
since  it  is  the  major  source  for  the  turbulent  kinetic  energy.  Comparing  Figs.  1  10  indi¬ 

cates  that  latent  heating  contributes  very  little  to  the  turbulent-kinetic-energy  generation 
in  this  particular  simulated  case,  which  is  stable  with  respect  to  cloud-top  entrainment  in¬ 
stability. 

h.  Partitioning  the  moiature  flux 

The  moisture  flux  and  its  partitioning  for  all  three  stratus  cases  are  given  in  Fig.  12. 
The  total  flux  should  be  linear  in  height  for  a  homogeneous,  quasi-steady  state.  We  as¬ 
sume  its  partitioned  fluxes  are  also  linear  in  height.  Since  radiation  and  latent  heating 
do  not  change  the  total  moisture  field,  the  partitioning  of  the  moisture  fiux  is  straightfor¬ 
ward.  The  ENT  curve  is  the  component  response  to  entrainment,  while  the  SFC  curve  is 
that  to  the  surface  moisture  flux.  The  sum  (dotted  lines)  of  these  two  curves  agrees  well 
with  the  wqT  profile  computed  from  the  LES  data  (solid  curves). 

6.  Summary  and  conclusions 

We  have  analyzed  four  large-eddy  simulations  of  idealized  nocturnal,  stratus-topped 
PBLs  and  the  clear  convective  PBL.  The  simulations  are  designed  so  we  can  approxi¬ 
mately  isolate  the  physical  processes  of  (1)  cloud-top  radiative  cooling,  (2)  entrainment, 
(3)  surface  heating,  and  (4)  latent  heating,  and  study  their  contributions  to  the  turbulent 
transports. 
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First,  u'e  showed  that  the  turbulent  transports  of  heat  and  moisture  within  the  mined 
stratus-topped  boundary  layers  are  well-expressed  by  the  thermodynamic  differences  be¬ 
tween  updrafts  and  downdrafts  of  turbulent  circulations.  This  is  the  basis  for  convective 
mass  flux  modeling,  which  has  recently  been  applied  to  the  buoyancy-driven  PBL.  We 
showed  that  the  convective  mass  flux  Me  =  Q.^<Tyj  is  a  good  approximation,  not  only  for 
clear  boundary  layers  as  suggested  by  Businger  and  Oncley  (1990),  but  also  for  the  heat 
and  moisture  transports  throughout  the  stratus-topped  mixed  layer. 

Therefore,  we  examined  the  contribution  of  each  physical  process  to  these  differences 
between  updrafts  and  downdrafts,  by  conditionid  sampling  of  the  large-eddy  simulated 
fields.  This  gives  us  qualitative  pictures  of  how  each  process  contributes  to  the  turbulent 
transports.  Based  on  that,  we  qusmtitatively  examined  the  linear  process  partitioning  of 
the  thermodynamic  energy  and  moisture  fluxes.  First  we  assumed  that  each  process  con¬ 
tributes  linearly  in  height  to  H  and  togr  fluxes.  We  then  showed  that  each  component  flux 
contributes  to  the  net  flux  consistently  (at  all  heights)  with  those  summarized  in  Section 
4c:  The  cloud-top  radiative  cooling  process  gives  a  positive  contribution  to  H,  entriun- 
ment  a  negative  one,  and  only  when  the  former  dominates  the  latter  can  the  .ff-flux  re¬ 
mains  positive  in  the  absence  of  surface  heating.  In  the  mid-PBL  of  the  ENT/RAD/SFC 
(dust  cloud)  case,  the  surface  heating  contributes  efficiently  to  more  than  half  of  the  H. 
The  surface  effect  contributes  about  twice  as  luge  to  the  net  moisture  flux  in  the  mid- 
PBL  of  the  ENT/RAD/SFC  case  comparing  to  the  ENT/RAD  case.  These  consistencies 
provide  a  justification  for  the  linear  process-peirtitioning  technique. 

We  also  showed  that  including  only  the  within-turbulent-layer  radiation  effect  (about 
90%  of  the  total  flux  jump)  does  the  sum  of  all  component  fluxes  give  an  H-Rux  profile 
that  agrees  well  with  that  computed  directly  &om  the  LES  data. 

Therefore,  the  S  profile  can  be  determined  from  three  quantities:  (1)  the  surface  flux, 
{whi)o,  (2)  the  entrainment  flux,  (tffhj),^,  and  (3)  the  bulJc-flux-jump  of  longwave  radiation 
across  the  whole  turbulent  layer,  AF,  which  excludes  the  about  10%  of  the  radiative  flux 
jump  that  lies  within  the  inversion  layer.  After  we  determine  H,  we  can  then  determine 
the  vertical  profile  of  whi  by  subtracting  J*Vatf  from  H,  with  Frad  m  a  prescribed  forcing. 
This  way,  we  include  the  effect  of  the  vertically-distributed  longwave  cooling  on  the  turbu- 
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lent  transports,  as  pointed  out  by  Randall  (1930b).  Lilly  and  Schubert  (1980),  and  Hanson 
(1987b). 

Similaxly,  we  can  determine  the  vertical  profile  of  the  total  moisture  flux  i37^  from 
Together,  whi  and  vjqr  determine  the  buoyancy  flux  profile. 

The  main  purpose  of  this  paper  is  to  show  methods  to  identify  the  contributions  of 
various  processes  to  the  updrafts  and  downdrafts,  and  to  the  flux  profiles,  given  the  entrai- 
ment  rate,  longwave  radiation  flux,  and  the  surface  flux.  Whether  the  process  partition¬ 
ing  technique  can  be  used  to  determine  the  entrainment  rate,  as  proposed  by  Stage  and 
Businger  (1981),  is  beyond  the  scope  of  this  paper  to  discuss. 

Our  discussions  here  are  for  a  highly-mixed  buoyancy-driven  PBL,  where  the  domi¬ 
nant  turbiilent  circulations  fill  the  whole  PBL.  When  the  stratus  cloud  layer  is  decoupled 
from  its  subcloud  layer  (because  of  internal  solar  heating  or  drizzling),  the  turbulent  circu¬ 
lations  no  longer  fill  the  whole  PBL  and  linear  partitioning  no  longer  applies. 

Also,  our  discussions  apply  only  to  stratus-topped  PBL  that  is  stable  with  respect  to 
the  cloud-top  entrainment  instability.  When  the  cloud-top  evaporative  cooling  effect  dom¬ 
inates,  there  is  likely  a  strong  entrainment /evaporation  interaction:  Stronger  entrainment 
leads  to  larger  cloud-top  evaporative  cooling  with  more  cooling  in  the  downdraft  branch. 
This  can  lead  to  stronger  turbulent  circulations  and  consequently  to  an  enhanced  entrain¬ 
ment.  Whether  this  positive  feedback  interaction  will  lead  to  a  breakup  of  the  stratus  deck 
is  stiU  in  debate  (e.g.,  Hanson,  1984;  Kuo  and  Schubert,  1988;  Siems  et  al.,  1990;  Albrecht, 
1991).  In  the  near  future,  we  will  analyze  the  large-eddy  simulation  of  a  stratus-topped 
PBL  that  is  unstable  with  respect  to  this  instability,  and  examine  this  feedback  mecha¬ 
nism. 
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FIGURE  CAPTIONS 


Figure  1  Vertical  distributions  of  the  buoyancy  fluxes  from  the  four  large-eddy  simula¬ 
tions  described  in  Table  1. 

Figure  2  Vertical  distributions  of  whi  (or  W3„  for  the  clear  case).  The  solid  curves  are  ob¬ 
tained  by  correlating  the  resolved-scale  vj  and  hi  (or  Sv)  of  the  large-eddy  simiilated 
fleld,  and  the  dashed  curves  are  from  equation  (2)  with  6=0.6. 

Figure  3  Same  as  Fig.  2,  but  for  w^. 

Figure  4  Vertical  distributions  of  the  mean  verticid  velocities  within  updrafts  (solid 

curves)  and  within  downdrafts  (dashed  curves),  deviating  from  the  “horizontal”  mean 
value. 

Figure  5  Same  as  Fig.  4,  but  for  the  virtual  temperature  fleld. 

Figure  6  Same  as  Fig.  4,  but  for  the  total  moisture  fleld. 

Figure  7  Same  as  Fig.  4,  but  for  the  liquid  water  fleld. 

Figure  8  Same  as  Fig.  4,  but  for  the  longwave  radiative  flux. 

Figure  9  Vertical  distributions  of  the  total  thermodynamic  energy  flux  H  (solid  curves) ^ 
the  longwave  radiative  flux  (dashed  curves),  and  the  hypothesized  longwave  radiation 
contribution  to  H  (dotted  ctirves). 

Figure  10  The  vertical  distributions  of  whi  from  the  large-eddy  simulations. 

Figure  11  Partitioning  of  the  total  thermodynamic  energy  flux  H  into  the  cloud-top  ra¬ 
diative  cooling  process  (RAD  curves),  the  entrainment  process  (ENT  curves),  and  the 
surface  heating  process  (SFC  curves).  Their  sum  is  shown  as  the  dotted  line. 

Figure  12  Partitioning  of  the  total  moisture  fliixes  wgx  into  the  entrainment  process 
(ENT  curves)  and  the  surface  moistening  process  (SFC  curves).  Their  sum  is  shown 
as  the  dotted  line. 


Table  1:  scaling  parameters 
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Figure  2  Vertical  distributions  of  w/ij  (or  ws„  for  the  clear  case).  The  solid  curves  are  ob 
tained  by  correlating  the  resolved-scale  w  and  hi  (or  s,)  of  the  large-eddy  simulated 
field,  and  the  dashed  curves  are  from  equation  (2)  with  6=0.6. 
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Figure  4  Vertical  distributions  of  the  mean  verticsd  velocities  within  updrsits  (sohd 

curves)  and  within  downdrafts  (dashed  curves),  deviating  from  the  “horizontal”  mesin 
value. 
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Figure  5  Same  as  Fig.  4,  but  for  the  virtual  temperature  field 
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Figure  6  Same  as  Fig.  4,  but  for  the  total  moisture  field 
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Figure  9  Vertical  distributions  of  the  total  energy  flux  H  (solid  curves),  the  longwave  ra¬ 
diative  flux  (dashed  curves),  and  the  hypothesized  longwave  radiation  contribution  to 
the  total  energy  flux  H  (dotted  curves). 
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Figure  10  The  rerticsl  distributions  of  whi  from  the  large-eddy  simulations. 
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